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Abstract 



The NA50 Collaboration has measured heavy-quarkonium production in p-A collisions at 
450 GeV incident energy (^/s = 29.1 GeV). We report here results on the production 
of the T states and of high- mass Drell-Yan muon pairs (m^^ > 6 GeV/c^). The cross- 
section at midrapidity and the A-dcpcndence of the measured yields are determined and 
compared with the results of other fixed-target experiments and with the available theoret- 
ical estimates. Finally, we also address some issues concerning the transverse momentum 
distributions of the measured dimuons. 
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1 Introduction 



The study of heavy quarkonium states is very important for our understanding of 
the physics of strong interactions and represents today one of the most challenging 
fields of application of QCD [1]. By studying quarkonium production in p-A col- 
lisions at fixed-target energies, various aspects of the theory of strong interactions 
can be addressed. In particular, such data can be used to constrain theoretical ap- 
proaches used for the calculation of the production cross-section in nucleon-nucleon 
interactions, such as the color-evaporation model (CEM) [2] or the nonrelativistic 
QCD formulation (NRQCD) [3]. Furthermore, by comparing the production data 
on various nuclear targets, the details of the color neutralization process of the qq 
pair can be investigated [4-6]. Finally, p-A data are an essential reference for the 
study of the anomalous suppression of quarkonia production predicted to occur in 
ultrarelativistic heavy- ion collisions [7-9]. 

To investigate these topics, the NA50 Collaboration performed a systematic study 
of quarkonium production on several nuclear targets with proton beams at 450 GeV/c 
incident momentum. Results on the charmonium states J/ip and detected through 
their decay into two muons, have already been published [10,11]. Thanks to the 
large integrated luminosity (~10 pb^^ for each target), the dimuon mass spectrum 
also shows a signal corresponding to the sum of the T(1S), T(2S) and T(3S) states 
(simply denoted as T in the following). 

Various results on T production in proton-induced collisions exist at fixed target 
and at ISR energies [12-25], but they are usually less precise than the corresponding 
charmonium results. Up to now, no data are available at 450 GeV {^/s = 29.1 GeV), 
and the A-dependence of the production cross-section has only been studied by the 
E772 experiment [22], at a higher incident energy. 

In this Letter we show results on the production of the T in proton-induced 
collisions on five nuclear targets (Be, Al, Cu, Ag, W). After applying the appropriate 
cuts, the total number of T meson events is of the order of 300. We compute, for the 
various targets, the cross-sections around midrapidity for the T, as well as for Drell- 
Yan muon pairs with invariant mass larger than 6 GeV/c^. We parametrize these 
cross-sections with the usual relation a a — ctnA'^, and determine the values ar and 
ckdy- We also determine the A-dependence of the ratio Bn^ar / cfdy-i a quantity less 
affected by systematical effects. Then, by combining the results for the different target 
nuclei, we obtain B ^^da / dycm\ycm=Oi ^ quantity that can be compared with the results 
of other experiments and with theoretical calculations. Finally, the dependences of 
(Pt) and (p^) on A and ^/s are discussed. 

2 Experimental setup and data analysis 

The setup of the NA50 experiment has been described in detail in various publica- 
tions [26]. It is based on a muon spectrometer composed of a 5 meter long hadron 
absorber, an air-core toroidal magnet and various sets of MWPCs and scintillator 
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hodoscopes. The luminosity measurement is achieved by means of three argon ioni- 
sation chambers placed along the beam path. They have been calibrated at low beam 
intensity and their linearity has been checked up to 10^° p/s [27]. The intensity of 
the 450 GeV proton beam dehvered by the SPS was about 3-10^ p/burst, with a 2.3 s 
spill. The thickness of the targets ranged from 30 to 50% of an interaction length. In 
these conditions, the dimuon trigger rate was ~10^/burst, with a dead time smaller 
than 5%. The efficiency etrig of the dimuon trigger system has been measured during 
data taking, using dedicated hardware, and ranges from 87 to 89% (see Table H)). 
About 50% of the triggered events contain a dimuon, which is reconstructed with 
the improved version of the reconstruction algorithm described in Refs. [11,28]. The 
dimuon reconstruction efficiency, e^^, is dominated by the efficiency of the MWPCs 
and is determined by the reconstruction program. Its precise value, used as input 
to obtain the cross-sections, depends on the data taking conditions and ranges from 
84 to 87% (see Table [T)). Various quality cuts are applied, to remove parasitic off- 
target events. Furthermore, in order to discard dimuons produced at the edge of the 
acceptance of the spectrometer, the kinematical domain of the analysis is restricted 
to —0.5 < Ucra < 0.5 in rapidity and to —0.5 < cos6'cs < 0.5 in the Collins-Soper 
polar angle of the muons [29]. We are finally left with about 1000 high mass events 
('^/i/i > 6 GeV/c^) for each target. 

Figure H shows the five invariant mass spectra, in the region m^^ > 6 GeV/c^. 
The T signal can be clearly seen on top of a continuum due to the Drell-Yan process. 
In this mass region the combinatorial background due to vr and K decays, usually 
calculated through the measured like-sign dimuon sample, turns out to be negligible. 
Muon pairs originating from Drell-Yan and from the T resonances are estimated by 
means of a fit to the invariant mass spectrum. The mass shapes for the various 
processes have been calculated through a Monte-Carlo simulation. The Drell-Yan 
events have been generated using a LO QCD calculation, with the GRV94LO set 
of parton distribution functions [30]; the distributions have been tuned directly 
on the data (see Section El for details) and, in the considered mass range, they do 
not significantly depend on the dimuon mass. The bottomonium states have been 
generated at their nominal PDG mass [31]. Due to the invariant mass resolution of 
the spectrometer, of the order of 0.4 GeV/c^ for dimuons with m^^ ~ 10 GeV/c^, 
the various T states cannot be resolved. We have, therefore, relied on previous 
measurements to fix the relative contributions of the three states visible in the dimuon 
channel. We have used the values T(2S)/T(1S) = 0.32 and T(3S)/T(1S) = 0.13, as 
measured in p-A collisions at 400 GeV [14], i.e. close to our energy. The results 
of Ref. [14] are anyway compatible with a high statistics measurement performed 
at higher energy [21]. For the bottomonium distribution, we have adopted the 
functional form 

r 21 6 

da/dpT oc pt/ [I + (pt/pq) \ , (1) 

used in Ref. [21]. The po parameter has been fitted on the measured pi distributions 
(see Section E)). 

For the rapidity, due to the narrow coverage of the set-up and to the relatively 
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small statistics, our data do not allow a precise determination of the T y-distribution. 
As a first guess, we have used in the acceptance calculation a gaussian distribution, 
centered at ?/cm=0, with ay = 0.37 ± 0.01. This value comes from a fit to the xp 
distribution for T production at 400 GeV (cTy = 0.35 ± 0.01) [14, 19]. The obtained 
value has then been logarithmically scaled to our incident energy. In Section El we 
will investigate in more detail the influence of the rapidity shape on our T results. 
Finally, the cos 9cs distributions for the T states have been generated according to the 
recent results of the E866 Collaboration [24], which found a negligible polarization for 
the T(1S) around midrapidity and a large transverse polarization for the T(2S) and 
T(3S) states. The generated Drell-Yan and T events have then been tracked through 
the apparatus and the accepted events have been reconstructed using the procedure 
applied to the experimental data. With the choice of the generation parameters 
described above, the acceptances for the various dimuon sources turn out to he Ady = 
21% (m^^ > 6 GeV/c2), At(is) = 25%, At(2S) = 25%, At(3S) = 26%. These 
acceptances refer to the phase space domain V defined by the cuts —0.5 < ?/cm < 0.5 
and —0.5 < cos 6*03 < 0.5. The error on the acceptances due to the uncertainty on 
the rapidity distribution is less than 1 %. 

Having calculated the shapes of the expected contributions, we fit each of the five 
invariant mass spectra with the function 

diV^- _ ,rDY^N^\,rr(dN^'''' ^^^ + 13 

dM dM \ dM ■ At(is) dM ' At(is) dM J 

(2) 

where N^"^ and A^^ are free parameters in the fit. To compensate for possible 
approximations in the description of the setup and in the mapping of the magnetic 
field we leave the position of the T(1S) peak as a further free parameter in the fit. The 
mass differences mx(2S) ~ ''^T(is) and mx(3s) — ''^T(iS') are anyway kept fixed to their 
PDG values [31]. We find that the shift from the nominal mass of the bottomonium 
states does not exceed 1% and has negligible consequences (< 1%) on the estimated 
T yield. 

The choice of the the fit region is dictated by the request of having a negligible 
contribution of the charmonia states. To satisfy this request, one could start the fit 
at mass values around 4.5 GeV/c^. However, due to the steep slope of the Drell-Yan 
mass distributions, the use of a large fit window results in a loss of weight of the low- 
statistics T region. In this way, small inaccuracies in the Monte-Carlo description 
of the Drell-Yan might influence the estimate of the T yield. Pushing down to 4.5 
GeV/c^ the starting point of the fit, the change in the calculated number of Drell-Yan 
and T events is only 10% on average, but the of the fits are worsened by about 
50%. Therefore, we have started our fits at m^^ = 6 GeV/c^. In this way, as can be 
seen in Fig. ^ the measured invariant mass spectra are well reproduced, with values 
of /ndf ranging from 0.9 to 1.5. The number of T mesons and the number of high 
mass Drell-Yan pairs are reported in Table ^ 

Finally, the Drell-Yan and T cross-sections, in the kinematical domain V defined 
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above, have been calculated through the formula 
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(3) 



where Ni is the number of detected events for each process £ is the luminosity, 
corrected for dead time, Ai is the acceptance, e^^ is the dimuon reconstruction effi- 
ciency and etrig is the efficiency of the dimuon trigger. For the T, we calculate the 
cross-section integrated over the T(1S), T(2S) and T(3S) states. The values of the 
various quantities relevant for the determination of the cross-section are summarized 
in Table CI 

3 Results 

In Fig. 121 we show, as a function of the mass number A, the Drell-Yan cross-sections 
divided by A. The points refer to the kinematical domain X>, and to the mass range 
m^^ > 6 GeV/c^. The error bars shown in Fig. El represent the quadratic combination 
of statistical and systematical errors. The systematical errors on the cross-section 
measurements are due to the errors on the determination of the incident proton flux 
and on the evaluation of etng and e^^ [11]. They range from 3.4 to 3.7% and are 
comparable to the statistical errors. The quantity aovIA corresponds to the cross- 
section per nucleon-nucleon collision. In absence of flnal state interactions, and if 
nuclear effects on the PDFs are negligible, a flat behaviour is expected. By fltting 
the points with the function 



we get aDy=0.98±0.02 {x^/ndf = 3.1), a value compatible with 1. Alternatively, 
imposing a=l results in the line shown in Fig. |21 with /"^df = 2.5. This result is 
in fair agreement with previous observations by NA50 in a lower mass range [10,11], 
showing that the Drell-Yan cross-section scales with the number of nucleon-nucleon 
collisions. 

Figure |21 shows the A-dependence of bottomonium production, through the ratio 
B^^ar /(Jdy, where B^^a-f refers to the sum of the IS, 2S and 3S states, weighted 
by their branching ratio into two muons. The quantity B^^a^ / ctoy is proportional 
to the bottomonium cross-section per nucleon-nucleon collision. Being the ratio of 
two measured cross-sections it is less affected by systematical errors. In the same 
figure we show the quantity B^^a^/A as a function of A. Again, the plotted errors 
represent a quadratic combination of statistical and systematical errors. The values 
of the various cross-sections and of their ratios is shown in Table |21 

The fits of the two sets of points to the usual A" parameterization give ar/DY = 
0.98 ±0.10 (xVnrf/=1.7) and ar = 0.98 ± 0.10 {x^ /ndf =1.7). Within the rather 
large statistical errors, our result indicates that the T is not strongly absorbed in 
the nuclear medium. In particular, by imposing a = 1 (solid lines in Fig. we 





(4) 
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can describe the data with x^/^df = 1.3 for T/DY and for ar- We find that our 
result is compatible with the only other direct determination of the A-dependence of 
T production, carried out by E772 [22] at 800 GeV incident energy. 

Assuming ar = 1, as suggested by our data, we averaged the results on the 
various targets and obtained for the T cross-section per nucleon-nucleon collision the 
value B^^ar = 0.27 ± 0.03 pb/nucleon. As a further check of possible systematic 
effects, we also obtained the same quantity through the analysis of a different data 
sample, taken by NA50 [10] at a beam intensity about one order of magnitude lower 
and containing 72±12 T events. We get B^^ar = 0.32±0.06 pb/nucleon. The two 
values of ux agree within errors, indicating that the high-luminosity data analyzed 
in this paper have been properly corrected for the various efficiency factors, whose 
evaluation is more delicate at high beam intensities. For completeness, we show 
in Fig. I3J as a function of A, B^^ar /ctdy and -B^j^crx for the low beam intensity 
sample. Performing a simultaneous fit of the data sets of Fig. |3] and |31 we get 
ay/oY = 0.98 ±0.09 {x'^ /ndf =0.9) and ax = 0.98 ± 0.08 {x^ /ndf =0.8). 

The T cross-section per nucleon-nucleon collision can be compared with previous 
measurements performed by various fixed-target experiments. However, past experi- 
ments usually quote the quantity Bf^^^dcr /dycm\ycm=o, and assumptions different from 
ours have sometimes been adopted in the cross-section calculations. In particular, 
past experiments evaluated their acceptances assuming an unpolarized production 
of the T states. Furthermore, in some cases they applied a correction to take into 
account the Fermi-motion of the target nucleons. In order to properly compare our 
results with the ones available in the literature, we have first re-calculated our accep- 
tances assuming no polarization for all the three bottomonium states. With respect 
to the results shown in Fig. El and 0] we get a 1.5% decrease in the value of the cross- 
section. The influence of Fermi motion has been taken into account by applying to 
our results (and to the ones available in the literature not yet corrected) the correc- 
tion proposed in Ref. [13]. At our center of mass energy, y/s = 29.1 GeV, the factor 
to apply to the results is 0.897. 

Concerning the A-dependence, all the other experiments also assumed a=l. 

A further element in the calculation of the midrapidity cross-section is the rapidity 
shape assumed for the T production. In fact, due to the limited y-coverage of the 
NA50 setup, our result may be sensitive to the chosen rapidity shape. We find 
that the rapidity distribution of the T mesons in our sample is compatible with a 
gaussian having > 0.30 (95% c.l.). Indeed, it is even compatible with a fiat y- 
distribution. The value for i?^^dcr/d?/cm|?/cm=o varies by ^^1"^% in these two extreme 
cases, and we consider this uncertainty as a contribution to the systematical error in 
the determination of the mid-rapidity cross- sect ion. 

The obtained value is Bfj_^da/dy\y^^=Q = 0.65 ± 0.05lHo 

pb/nucleon, where the 

quoted asymmetric error is due to the uncertainty on the T rapidity distribution as 
derived from our data. In Fig.^lwe present the i/i-dependence of the T cross-section 
per nucleon-nucleon collision. The NA50 point follows the general trend defined by 
previously available results. We also plot in Fig. El the results of a NLO calculation 
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based on the Color Evaporation Model (CEM) [33]. The curves are obtained by fit- 
ting the data with various combinations of the CEM parameters, namely the mass of 
the b-quark, the renormalization/factorization scale and the choice of the PDF set. It 
can be seen that the ^/s dependence of the T cross-section is well reproduced. For the 
other model commonly used in the study of heavy quarkonia cross-sections, NRQCD, 
no up-to-date calculation exists in the literature. The only available prediction [34] 
gives a value of the order of 50 pb/nucleon for the total T cross-section at y/s = 29.1 
GeV. By integrating our differential cross-section over y and cosOcs^ and correcting 
for the branching ratios, we get cr^^ = 29±2l4^ pb/nucleon. The asymmetric system- 
atic error is dominated by the uncertainty in the rapidity distribution derived from 
our data. If we rely on the data from previous experiments for the determination of 
(tJ the systematic error is reduced to ±2 pb/nucleon. 

Finally, we have investigated the A-dependence of the dimuon transverse momen- 
tum distributions. In Figs. IHl and d we show, for the various systems, the px spectra 
observed in the mass intervals 4.5 < m^^ < 8 GeV/c^ and 8.6 < m^^ < 11.6 
GeV/c^, respectively. In Fig. IHl the Drell-Yan distributions, generated according 
to Eq. Hand filtered through the Monte-Carlo simulation of the detector, have been 
fitted to the data with po as a free parameter. In Fig. E|the spectra have been fitted 
with a superposition of Drell-Yan and T dimuons. For Drell-Yan (dashed lines) the 
normalizations have been obtained from the fit to the mass distributions of Fig. Q 
while the shapes have been extrapolated from the pt spectra measured in the lower 
mass range. The T events have been generated following Eq. ^ and po has been used 
as a free parameter in the fits. As a result, we find that for the two mass intervals 
under study our fits nicely reproduce the observed distributions. In Table El we 
summarize the po values obtained for the various systems, for both Drell-Yan and T. 
The calculated (px) and (px) are also shown. To test the stability of our results on 
the choice of the assumed Drell-Yan px distributions, we have performed a similar 
study, fitting our data to a PYTHIA calculation, performed with the GRV94LO set 
of parton distribution functions. The value of the intrinsic {k±) of partons inside 
hadrons, used in PYTHIA, has a strong infiuence on the Drell-Yan px distribution. 
It has been tuned directly on the data, obtaining values ranging from 0.88 to 0.90 
GeV/c. We find that the shape of the PYTHIA px distribution still reproduces our 
data, with a slightly worse x^/^df. The calculated (px) differ by less than 3% (7% 
for (px)) from the result obtained using Eq. [T] 

According to [35,36], for heavy quarkonia and Drell-Yan production, initial state 
parton scattering is responsible for the increase of (px) observed in p-A collisions 
with respect to pp. Such an increase is approximately linear as a function of L, 
the mean length of nuclear matter crossed by the incoming parton, and its size is 
smaller (by a factor 4/9) if the scattering process is initiated by a quark rather than 
by a gluon. In Ref. [32] it was shown that the transverse momentum distribution of 
the J/ip mesons produced in p-A and A-A collisions can be fitted with the function 
(p2^ = (p2,)pp + ■ L, with ttgN = 0.077±0.002 {GeV/c)Hm-\ In Fig. H we show 
the L-dependence of (p^) for Drell-Yan dimuons in the mass region 4.5 < m^^ < 8 
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GeV/c^, where the initial hard parton is a quark rather than a gluon. A hnear increase 
of {p\) with L is visible. By using the same kind of fit, we get agAr=0.021±0.019 
(GeV/ c)^fm~^, a value significantly lower than QgN and, within the rather large errors, 
in agreement with the expected 4/9 factor. For T production, where both qq and 
gg initial states contribute [34], we would expect for the (p^) increase vs. L an 
intermediate slope between the ones observed for Drell-Yan and J/'?/'. Unfortunately, 
the available statistics (see Fig. |HI) does not allow to draw quantitative conclusions. 
Finally, in Fig. El we compare, as a function of -^i, the measured (px) (averaged over 
the various target nuclei) with the results of other experiments. For Drell-Yan, a 
clear increase of (px) with y/s is observed, while the effect is less important for the 
T. 

4 Conclusions 

The NA50 experiment measured high-mass {m^^ >6 GeV/c^) Drell-Yan dimuons 
and, for the first time, T production in p-A collisions at = 29.1 GeV. By fitting the 
cross-section results with the usual yl" parameterization, we get aov = 0.98 ± 0.02 
and ax = 0.98 ± 0.08. Within the rather large experimental errors, this result indi- 
cates the absence of strong absorption effects for the bottomonium states in nuclear 
matter. The measured cross-section at mid-rapidity is in good agreement with the 
results of other experiments in the fixed target energy range and with the available 
theoretical calculations. The A-dependence of {p\) indicates, for Drell-Yan produc- 
tion, an effect compatible with quark scattering in the initial state. 
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Nt 


Ndy (nv^ > 6 GeV/c2) 


C (pb-i) 




^trig 


p-Be 


26±9 


485±23 


52.2 


0.86 


0.87 


p-Al 


82±13 


901±31 


30.2 


0.87 


0.88 


p-Cu 


67±13 


1110±34 


17.1 


0.84 


0.89 


p-Ag 


65±14 


1243±37 


9.9 


0.86 


0.89 


p-W 


63±12 


820±30 


4.7 


0.84 


0.87 



Table 1: Quantities used in the cross-section calculation. 





o-^y/A (pb) 


B^^ar/k (pb) 


B f_if_,crr / a DY 


p-Be 


5.01±0.24±0.18 


0.218±0.072±0.008 


0.044±0.015 


p-Al 


5.36±0.19±0.25 


0.399±0.063±0.019 


0.075±0.012 


p-Cu 


4.95±0.15±0.19 


0.245±0.046±0.009 


0.049±0.010 


p-Ag 


5.61±0.17±0.19 


0.239±0.052±0.008 


0.043±0.009 


p-W 


4.57±0.17±0.17 


0.283±0.055±0.010 


0.062±0.013 



Table 2: Drell-Yan (for the mass region m^^ >6 GeV/c^) and bottomonium cross- 
sections, and their ratio. The first error is statistical, the second systematical. The 
quoted systematical errors cancel out in the ratio B^^ay /acY- 





Po (DY) 
(GeV/c) 


Po (T) 
(GeV/c) 


(Pt) (DY) 
(GeV/c) 


(Pt) (T) 
(GeV/c) 


(pD (DY) 
(GeV/c)2 


{pD (t) 

(GeV/c)2 


p-Be 


2.77±0.04 


3.4±0.6 


1.19±0.03 


1.5±0.2 


1.91±0.05 


2.9±1.0 


p-Al 


2.75±0.03 


3.3±0.4 


1.18±0.01 


1.4±0.2 


1.89±0.04 


2.7±0.6 


p-Cu 


2.80±0.02 


3.0±0.4 


1.20±0.01 


1.3±0.2 


1.96±0.04 


2.3±0.7 


p-Ag 


2.78±0.02 


2.6±0.3 


1.19±0.01 


l.liO.l 


1.93±0.03 


1.7±0.4 


p-W 


2.80±0.03 


3.0±0.4 


1.20±0.01 


1.3±0.2 


1.96±0.04 


2.3±0.7 



Table 3: po (see Eq. H}, (pt) and for Drell-Yan (4.5 < m^^ < 8.0 GeV/c^) and 
T. The quoted errors for Drell-Yan are purely statistical, while the values for the T 
include a systematical error due to the uncertainty in the extrapolation of the Drell- 
Yan yield into the T region. The total error is anyway dominated by the statistical 
contribution, due to the low T statistics. 
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M (GeV/c^) M (GeV/c^) 

Figure 1: The five p-A opposite-sign muon pair mass spectra, in the mass range 
m^^ >6 GeV/c^. The sohd hne is the resuh of the fit with the function described 
in Eq. |21 The dashed hne represents the DreU-Yan process, the dotted hues the 
contributions of the various bottomonium states, the dot-dashed hne the sum of the 
three T states. 
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Figure 2: The Drell-Yan cross-section, relative to the mass region m^^ >6 GeV/c^, 
divided by the mass number A. The hne represents the result of a fit to the points, 
according to the function = a^y " ^■ 
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Figure 3: The ratio B^^ar /(Jdy (closed circles) and the cross-section per 
nucleonS^^crr/^ (open circles), as a function of A. 
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Figure 4: The ratio B^^ar /o-dy (closed squares) and tiie cross-section per nucleon 
B^^ar/A (open squares), as a function of A, for the low- intensity data sample. 
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Figure 5: T cross-section at midrapidity as a function of y^. The lines represent the 
results of NLO CEM calculations [33]. The solid, dashed and dot-dashed lines are 
obtained with the MRST HO PDF set with, respectively, nib = /i = 4.75 GeV/c^, 
nib — A*/2 = 4.5 GeV/c^, and nib — 2/1 — 5 GeV/c^. The dotted hne is obtained 
with the GRV HO PDF set with = // = 4.75 GeV/c^ It is assumed that 
IJ' — IJ'R — IJ'F, where /ir is the renormalization scale and /ip is the factorization scale. 
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Figure 6: The measured p-A opposite-sign dimuon px spectra in the invariant mass 
region 4.5 < m^^ < 8.0 GeV/c^. The hues represent the best fit to the data obtained 
with the functional form shown in Eq. ^ 
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Figure 7: The measured p-A opposite-sign dimuon px spectra in the invariant mass 
region 8.6 < m^^ < 11.6 GeV/c^. The dashed hues represent the contribution of the 
Drell-Yan process, the sohd hues the sum of the T and Drell-Yan yields (see text for 
details). 
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Figure 8: The L-dependence of the average transverse momentum squared for Drell- 
Yan dimuons in the mass range 4.5 < m^^ < 8 GeV/c^ (left plot) and for T (right 
plot). For Drell-Yan, the result of a linear fit is also shown. 
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Figure 9: The average transverse momentum for Drell-Yan (a) and T (b) produc- 
tion as a function of y/s. For Drell-Yan, the points correspond to m ~ 0.22 • y/s 
(corresponding to = m/y^ ~ 0.22). 
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